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Abstract 

Searches for heavy Majorana neutrinos in B~ decays in final states containing 
hadrons plus a fi~ pair have been performed using 0.41 fb _1 of data collected 
with the LHCb detector in proton-proton collisions at a center-of-mass energy of 
7 TeV. The D + n~ fi~ and D* + \T~ fj," final states can arise from the presence of 
virtual Majorana neutrinos of any mass. Other final states containing ir + , Df, or 
* D°ir + can be mediated by an on-shell Majorana neutrino. No signals are found and 

upper limits are set on Majorana neutrino production as a function of mass, and 
also on the B~ decay branching fractions. 
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1 Introduction 



Leptons constitute a crucially important sector of elementary particles. Half of the leptons 
are neutrinos. Yet we do not know if they are Dirac or Majorana particles, the latter 
case characterized by being their own antiparticles pp. Since the observation of neutrino 
oscillations has indisputably established that neutrinos have non-zero mass, it is possible 
to distinguish the two types experimentally. Finding neutrinoless double /3 decay has long 
been advocated as a premier demonstration of the possible Majorana nature of neutrinos 
[2] . The Feynman diagram is shown in Fig. [I] We also show the fundamental quark and 
lepton level process. An impressive lower limit from neutrinoless double decays in nuclei 
has already been obtained on the half-life of C(10 25 ) years [3] for coupling to e~. 




Figure 1: (a) Diagram of neutrinoless double /3 decay when two neutrons in a nucleus 
decay simultaneously, (b) The fundamental diagram for changing lepton number by two 
units. 

Similar processes can occur in B~ decays. The diagram is shown in Fig. |2^a). In 
this reaction there is no restriction on the mass of the Majorana neutrino as it acts as 
a virtual particle. In this paper, unlike in neutrinoless double beta decays, a like-sign 
dimuon is considered rather than two electrons. The only existing limit is from a recent 
Belle measurement jl] using the B~ — > D + fi~fi~ channel. We consider only final states 
where the cd pair forms a final-state meson, either a D + or a D* + , so the processes we 
are looking for are B~ — > n~ ji~ . In this paper mention of a specific reaction also 
implies inclusion of the charge conjugate reaction. 

There are other processes involving 6-quark decays that produce a light neutrino that 
can mix with a heavy neutrino, designated as N . The heavy neutrino can decay as 
N — > W + pT . In Fig. |2](b) we show the annihilation processes B~ — > 7r + (.D+)//~/i~, where 
the virtual W + materializes either as a tt + or Df. These decays have been discussed in 
the literature [5J E] . 

We note that it is also possible for the B~ — > £)(*)+ pr fjr decay modes shown in 
Fig. [2]^a) to proceed by a Cabibbo suppressed version of the process in Fig. [2^b) where 
the virtual W + forms D^ + . Similarly, the decay modes shown in Fig. unb) could be 
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Figure 2: Feynman diagrams for B decays involving an intermediate heavy neutrino (JV). 
(a) B- -> D^+fi-fi-, (b) 5- -> 7r+(D+)//-/i-, and (c) B~ ->■ D°n + fi~ fi~ . 

produced via Cabibbo suppressed versions of the process in Fig. [2(a). Here the 

final state requires a 6 — > u quark transition while for the _D+/i _ /i~ final state, one of the 

virtual must couple to as quark rather than a d. 

The lifetimes of N are not predicted. We assume here that they are long enough that 
the natural decay width is narrower than our mass resolution which varies between 2 
and 15 Me depending on mass and decay mode. For B~ — > Tr + /i~/j l '~, we can access 
the Majorana mass region between approximately 260 and 5000 MeV while for B~ — > 
Dfn~fi~, the Majorana mass region is between 2100 and 5150 MeV. In the higher mass 
region, the W + may be more likely to form a meson than a 7r + . The B~ — > tt + {i~h~ 
search was first performed by Mark-II j7j and then by CLEO [HJ. LHCb also performed a 
similar search using a smaller 0.04 fb _1 data sample j9| giving an upper limit of 5.8 x 10~ 8 
at 95% confidence level (CL). The decay of B~ — > D+yT \i~ has never been investigated. 

Finally, in Fig. |2](c) we show how prolific semileptonic decays of the B~ can result 
in the D°tt + [i~ final state. This process has never been probed [TU]. We benefit from 
the higher value of the CKM coupling \V c b\ relative to |V^,| in the annihilation processes 
shown in Fig. [2(b). The accessible region for Majorana neutrino mass is between 260 
and 3300 MeV. For all the modes considered in this paper, we search only for decays 
with muons in the final state, though electrons, and r leptons in cases where sufficient 
energy is available, could also be produced. Searches have also been carried out looking 
for like-sign dileptons in hadron collider experiments |llj . 

1 In this paper we use units where the speed of light, c, is set equal to one. 
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2 Data sample and signal 



We use a data sample of 0.37 fb -1 collected with the LHCb detector p2] in the first half 
of 2011 and an additional 0.04 fb _1 collected in 2010 at a center-of-mass energy of 7 TeV. 

The detector elements are placed along the beam line of the LHC starting with the 
vertex detector, a silicon strip device that surrounds the proton-proton interaction region 
having its first active layer positioned 8 mm from the beam during collisions. It provides 
precise locations for primary pp interaction vertices, the locations of decays of long-lived 
particles, and contributes to the measurement of track momenta. Further downstream, 
other devices used to measure track momenta include a large area silicon strip detector 
located in front of a 4 Tm dipole magnet, and a combination of silicon strip detectors 
and straw-tube drift chambers placed behind. Two Ring Imaging Cherenkov (RICH) 
detectors are used to identify charged hadrons. An electromagnetic calorimeter is used 
for photon detection and electron identification, followed by a hadron calorimeter, and a 
system that distinguishes muons from hadrons. The calorimeters and the muon system 
provide first-level hardware triggering, which is then followed by a software high level 
trigger. 

Muons are triggered on at the hardware level using their penetration through iron and 
detection in a series of tracking chambers. Projecting these tracks through the magnet 
to the primary event vertex allows a determination of their transverse momentum, pp- 
Events from the 2011 data used in this analysis were triggered on the basis of a single 
muon having a greater than 1480 MeV, or two muons with their product greater 
than 1.69 GeV 2 . To satisfy the higher level trigger, the muon candidates must also be 
detached from the primary vertex. 

Candidate B~ decays are found using tracking information, and particle identification 
information from the RICH and muon systems. The identification of pions, kaons and 
muons is based on combining the information from the two RICH detectors, the calorime- 
ters and the muon system. The RICH detectors measure the angles of emitted Cherenkov 
radiation with respect to each charged track. For a given momentum particle this angle is 
known, so a likelihood for each hypothesis is computed. Muon likelihoods are computed 
based on track hits in each of the sequential muon chambers. In this analysis we do 
not reject candidates based on sharing hits with other tracks. This eliminates a possible 
bias that was present in our previous analysis [2]. Selection criteria are applied on the 
difference of the logarithm of the likelihood between two hypotheses. The efficiencies and 
the mis-identification rates are obtained from data using K$, D* + — > tt + D°, D° — > K~n + 
and J/ip — > event samples that provide almost pure pion, kaon, and muon sources. 

Efficiencies and and rejection rates depend on the momentum of the final state par- 
ticles. For the RICH detector generally the pion or kaon efficiencies exceed 90% and the 
rejection rates are of the order of 5% [13]. The muon system provides efficiencies exceeding 
98% with rejection rates on hadrons of better than 99%, depending on selection criteria 
[T4] . Tracks of good quality are selected for further analysis. In order to ensure that tracks 
have good vertex resolution we insist that they all have p^ > 300 MeV. For muons this 
requirement varies from 650— 800 MeV depending on the final state. All tracks must be 



3 



Table 1: Charm and charmonium branching fractions 



Particle Final state Branching fraction (%) 



D° 


K~7T + 


3.89±0.05 [3] 


D+ 




9.14±0.20 [3] 


Dt 




5.50±0.27 p3] 


D *+ 


1T + D° 


67.7±0.5 [3] 


if;(2S) 


7T + 7T~J/lf) 


32.6±0.5 


j/i/j 




5.93±0.06 j3] 



inconsistent with having been produced at the primary vertex closest to the candidate B~ 
meson's decay point. The impact parameter (IP) is the minimum distance of approach 
of the track with respect to the primary vertex. Thus we form the IP \ 2 by testing the 
hypothesis that the IP is equal to zero, and require it to be large; the values depend on 
the decay mode and range from 4 to 35. 

3 Normalization channels 

Values for branching fractions will be normalized to well measured channels that have the 
same number of muons in the final state and equal track multiplicities. The first such 
channel is B~ -> J/ipK~. Its branching fraction is B(B~ -> J/ipK~) = (1.014 ±0.034) x 
10~ 3 [3J. We use the J /if) — > fi + fi~ decay mode. The product branching fraction of this 
normalization channel is (6.013 ± 0.021) x 10~ 5 , and is known to an accuracy of ±2%. 
The charm meson decay modes used in this paper are listed in Table [TJ along with their 
branching fractions and those of the charmonium decays in the normalization channels. 

To select the J/if)K~ normalization channel, the pi requirement is increased to 
1100 MeV for the K~ and 750 MeV for the muons. To select B~~ candidates we further 
require that the three tracks form a vertex with a \ 2 < 7, and that this B~ candidate 
points to the primary vertex at an angle not different from its momentum direction by 
more than 4.47 mrad, and that the impact parameter \ 2 of the B~ is less than 12. The 
same requirements will be used for the selection. The total efficiency for fi + {i~K~ 

is (0.99±0.01)%, where the fi + come from J/ip decay. 

The invariant mass of K~fi + fi~ candidates is shown in Fig. [3^a). In this analysis the 
invariant mass is required to be within 50 MeV of the J /if) mass. We use a Crystal 
Ball function (CB) to describe the signal [16J, a Gaussian distribution for the partially 
reconstructed background events, and a linear distribution for combinatorial background. 
The CB function provides a convenient way to describe the shape of the distribution, 
especially in the mass region below the peak where radiative effects often produce an 
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excess of events that falls away gradually, a so called "radiative tail" . The CB function is 



f(m;a,n,m ,cr) 



exp 



(m— mo)- 

m— mn 
a 



for > -a 

(J 

for < -a 

a — 



(1) 



where 



.4 



• exp 



a 



n 
\ot\ 



a . 



The measured mass of each candidate is indicated as m, while mo and a are the fitted 
peak value and resolution, and n and a are parameters used to model the radiative tail. 
We use the notation a in the rest of this paper to denote resolution values found from CB 
fits. 

Using an unbinned log-likelihood fit yields 47,224±222 B~ — > J j ipK~ events. Within 
a ±2o" signal window about the peak mass, taken as the signal region, there are 44,283 
of these events. The number of signal events in this window is also determined using the 
total number of events and subtracting the number given by the background fit. The 
difference is 119 events, and this is taken as the systematic uncertainty of 0.3%. The 
width of the signal peak is found to be 19.1 ±0.1 MeV. Monte Carlo simulations are based 
on event generation using Pythia [17], followed by a Geant-4 [18] based simulation of 
the LHCb detector [19J. The J/ipK~ mass resolution is 20% larger than that given by 
the LHCb simulation. All simulated mass resolutions in this paper are increased by this 
factor. 
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Figure 3: Invariant mass of (a) candidate J /ipK~ decays, and (b) candidate J/ipK~7i + 'K~ 
decays. The data are shown as the points with error bars. Both the partially reconstructed 
background and the combinatorial background are shown, although the combinatorial 
background is small and barely visible. The solid curve shows the total. In both cases 
the candidate fi + fi~ is required to be within ±50 MeV of the J/ip mass, and in (b) the 
dimuon pair is constrained to have the J/ip mass. 
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For final states with five tracks, we change the normalization channel to B — > 
ip(2S)K~, with ip(2S) — > 7r + 7r~ J/if>, and J/if> — > /i + /U~. The branching fraction for 
this channel is B(B~ -> ifj(2S)K~) = (6.48 ± 0.35) x 1CT 4 [3J. Events are selected using 
a similar procedure as for J/ifjK~ but adding a 7r + 7r - pair, that must have an invariant 
mass when combined with the J /if} which is compatible with the if>(2S) mass, and that 
forms a consistent vertex with the other B~ decay candidate tracks. The total efficiency 
for n + ji~ it + it~ K~ is (0.078±0.002)%, without inclusion of the if)(2S) or J/if> branch- 
ing fractions. The B~ candidate mass plot is shown in Fig. |3]^b). Here the pair 
is constrained to the J/if> mass. (In what follows, whenever the final state contains a 
ground-state charm meson, its decay products are constrained to their respective charm 
masses.) 

The data are fitted with a CB function for signal, a Gaussian distribution for partially 
reconstructed background and a linear function for combinatorial background. There are 
767±29 signal events in a ±2cr window about the peak mass. The difference between 
this value and a count of the number of events in the signal region after subtracting the 
background implies a 0.7% systematic uncertainty on the yield. 



4 Analysis of B D + \x \i and /x 

Decay diagrams for B~ — > n~ n~ are shown in Fig. |2](a). Since the neutrinos are 
virtual, the process can proceed for any value of neutrino mass. It is also possible for these 
decays to occur via a Cabibbo suppressed process similar to the ones shown in Fig. (2Tb), 
where the virtual W + materializes as a cd pair. If this occurred we would expect the 
Cabibbo allowed Z)+ \T [T final state to be about an order of magnitude larger. The search 
for Majorana neutrinos in this channel are discussed in Section [6j The D + — > K~tt + it + 
and D* + — > tt + D°, D° — > K~tt + channels are used. The decay products of the D + and 
D° candidates are required to have invariant masses within ±25 MeV of the charm meson 
mass, and for D* + candidate selection the mass difference m(ir + K~ 7r + ) — m(K~7T + ) is 
required to be within ±3 MeV of the known D* + — D° mass difference. 

The fi~ candidate mass spectra are shown in Fig. [Zj No signals are apparent. 

The B~ mass resolution is 15.7±0.5 MeV for the D + channel and 14.1±0.6 MeV for the 
D* + channel. The background has two components, one from mis-reconstructed B decays 
that tends to peak close to the B~ mass, called "peaking backgrounds" , and random track 
combinations that are parameterized by a linear function. To predict the combinatorial 
background in the signal region we fit the data in the sidebands with a straight line. 
In the D + mode we observe six events in the signal region, while there are five in the 
D* + mode. The combinatorial background estimates are 6.9±1.1 and 5.9±1.0 events, 
respectively. Peaking backgrounds are estimated from mis-identification probabilities, 
determined from data, coupled with Monte Carlo simulation. For these two channels 
peaking backgrounds are very small. The largest, due to B~ — > D + tt~tt~, is only 0.04 
events. 

The total efficiencies for D + /i' pT and D* + fx~fx~ are (0.099 ± 0.007)% and (0.066 ± 
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0.005)%, respectively; here the charm branching fractions are not included. The system- 
atic errors are listed in Table [2] for this mode and other modes containing charm mesons 
that will be discussed subsequently. Trigger efficiency uncertainties are evaluated from 
differences in the 2010 and 2011 data samples. The largest systematic uncertainties are 
due to the branching fractions of the normalization channels and the trigger efficiencies. 
The uncertainty on the background is taken into account directly when calculating the 
upper limits as explained below. Other uncertainties arise from errors on the charmed 
meson branching fractions. For these final states the uncertainty due to different final 
state track momenta with respect to the normalization mode are very small, on the order 
of 0.2%. Other channels have uncertainties due to varying efficiencies as a function of 
Majorana mass, and these are entered in the row labeled "Efficiency modeling". The 
detector efficiency modeling takes into account the different acceptances that could be 
caused by having different track momentum spectra. For example, the track momenta 
depend on the Majorana neutrino mass for on-shell neutrinos. These uncertainties are 
ascertained by simulating the detector response at fixed Majorana masses and finding the 
average excursion from a simple fit to the response and the individually simulated mass 
points. This same method is used for other modes. 
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Figure 4: Invariant mass spectrum for (a) B — > D + fi fi candidates, and (b) B — > 
D* + fi~fi~ candidates. The solid lines show the linear fits to the data in the mass sidebands. 



To set upper limits on the branching fraction the number of events iV bs within ±2<r 
of the B~ mass are counted. The distributions of the number of events (N) are Poisson 
with the mean value of (S + B), where S indicates the expectation value of signal and B 
background. For a given number of observed events in the signal region, the upper limit 
is calculated using the probability for iV < N^: 



P{N < N< 



obs J 



(S + B 



N<N n 



(2) 



A limit at 95% CL for branching fraction calculations is set by having P(N < N t, s ) = 0.05. 
The systematic errors are taken into account by varying the calculated S and B, assuming 
Gaussian distributions. 
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Table 2: Systematic uncertainties for B — > DXji ji modes. 



Source Systematic uncertainty (%) 

Common to all modes 



B(B ip(2S)K ) 

LjytpyZ,D J — r J 1 IfJTl 7T J 






5.4 
1 .o 
1.0 




Uncertainty in signal shape 






3.0 




Yield of reference channel 






0.7 




PL PID 






0.6 




Mode specific 


Dt 


_D°7T 4 


D + 


D*+ 


Trigger 


4.9 


9.3 


5.5 


4.8 


Efficiency modeling 


10.0 


6.7 






PID (K/n) 


1.0 








Charm decay £>'s 


4.9 


1.3 


2.2 


1.5 


Total 


13.8 


13.2 


8.8 


8.2 



The upper limits on the branching fractions at 95% CL are measured to be 

B{B~ D+ft-fj,-) < 6.9 x 10" 7 and 
B(B- D* + fi-fi-) < 2.4 x 10" 6 . 

The limit on the D + channel is more stringent than a previous limit from Belle of 1 x 10~ 6 
at 90% CL and the limit on the D* + channel is the first such result. 

5 Analysis of B — >- 7r + /i /i 

The selection of events uses the same criteria as described for J / ipK~ in Section|3j 

except for like-sign rather than opposite-sign dimuon charges and pion rather than kaon 
identification. The invariant mass distribution of candidates is shown in Fig. |5j 

The mass resolution for this final state is 20.3±0.2 MeV. An interval of ±2cr centered on 
the B~ mass is taken as the signal region. There are 7 events in the signal region, but no 
signal above background is apparent. The peaking background, estimated as 2.5 events, is 
due to misidentified B~ — > J /ipK~ or J/ip7i~ decays; the shape is taken from simulation. 
The combinatorial background is determined to be 5.3 events from a fit to the n + fi~fi~ 
mass distribution excluding the signal region. The total background in the signal region 
then is 7.8 ± 1.3 events. 

Since the putative neutrinos considered here decay into 7r + /i - , and are assumed to 
have very narrow widths, more sensitivity is obtained by examining this mass distribution, 
shown in Fig. [6| for events in the B~ signal region. There is no statistically significant 
signal at any mass. There are three combinations in one mass bin near 2530 MeV; 
however two of the combinations come from one event, while it is possible to only have one 



8 



> - 

<D " 

§4- 
in - 

tn '. 
c 3- 

<D - 

> - 

LU 

2 n 



1 - 



T 1 1 1 1 1 1 1- 

LHCb 



5100 



- I 1 1 1 1 — 

Peaking Background 
Combinatorial Background 



5200 5300 5400 

Invariant mass of 7t + |i M" (MeV) 



Figure 5: Invariant mass distribution of 7r + /i \l . The estimated backgrounds are also 
shown. The curve is the sum of the peaking background and the combinatoric background. 



Majorana neutrino per B~ decay. Upper limits at 95% confidence level on the existence 
of a massive Majorana neutrino are set at each 7r + /i~ mass by searching a signal region 
whose width is ±3<7tv, where <jn is the mass resolution, at each possible Majorana neutrino 
mass, M^. This is done in very small steps in 7r + /i~ mass and so produces a continuous 
curve. If a mass combination is found anywhere in the ±3<tat interval it is considered 
as part of the observed yield. To set upper limits the mass resolution and the detection 
efficiency as a function of 7r + /i~ mass need to be known. Monte Carlo simulation of the 
mass resolution as a function of the Majorana neutrino mass is shown in Fig. [7j along 
with resolutions of other channels. The overall efficiencies for different values of Mn are 
shown in Fig. [8j A linear interpolation is used to obtain values between the simulated 
points. 

Many systematic errors in the signal yield cancel in the ratio to the normalization 
channel. The remaining systematic uncertainties are listed in Table [3j The largest sources 
of error are the modeling of the detector efficiency (5.3%) and the measured branching 
fractions B (B~ -> J/ipK~) (3.4%), and B (J/ip -> (1-0%). 

To set upper limits on the branching fraction, the number of events iV b s at each 
value (within ±3ctjv) are counted, and the procedure described in the last section applied. 
Estimated background levels are taken from Fig. [6j Figure p)[a) shows the upper limit on 
B(B~ — > 7i + fi~fi~) as a function of Mn at 95% CL. For most of the neutrino mass region, 
the limits on the branching ratio are < 8 x 1CT 9 . Assuming a phase space decay of the 
B~ we also determine 

B(B- -»■ 7T+/T/0 < 1.3 x KT 8 at 95% CL. 
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Figure 6: Invariant mass distribution of 7r + /i~ in the ±2cr region of the B~ mass with 
both peaking and combinatorial background superimposed. The peaking background at 
3100 MeV is due to misidentified B~ — > J/ipX decays. There are two combinations per 
event. 




Majorana neutrino mass (MeV) 

Figure 7: Majorana mass resolutions for the three B~ decays as a function of Majorana 
mass. 
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Figure 8: Detection efficiencies for the three B decays as a function of Majorana mass. 
Charm meson decay branching fractions are not included. 



Table 3: Systematic uncertainties for B — > 7i + /j fi measurement. 



Selection criteria 


Systematic uncertainties (%) 


K/ir PID 


1.0 


fi PID 


0.6 


Muon selection 


0.6 


Trigger 


1.0 


Yields of reference channel 


0.4 


Efficiency modeling 


5.3 


B(B~ -> J/ipK-) 


3.4 


B(J/ip -> 


1.0 


Total 


6.7 



These limits improve on previous results from CLEO (< 1.4 x 10 6 at 90% CL) [8], and 
LHCb (< 5.8 x 10- 8 at 95% CL) ^. 

6 Analysis of B~ 

The process B~ — > D^jjTyT is similar to B~ — > 7T + fi~fi~, with the difference being 
that the heavy neutrino can decay into D^jjT. Here we consider only — > K + K~n + 
decays. Our analysis follows a similar procedure used for the tt + [i~{i~ channel. Candidate 
— > K + K~-k + decays are selected by having an invariant mass within ±25 MeV of the 
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Figure 9: Upper limits at 95% CL as a function of the putative Majorana neutrino mass, 
(a) for B(B~ — > n + n~ n~) as a function of the ir + fi~ mass, (b) for B{B~ — > D+fi~fi~) as 
a function of the Df[i~ mass, and (c) for B{B~ — > D°tt + fi~ fi~) as a function of the ir + fi~ 
mass. 



Df mass. A Majorana neutrino candidate decay is then looked for by by having the 
candidate decay tracks form a vertex with an opposite-sign muon candidate. Then 
this neutrino candidate must form a vertex with another muon of like-sign to the first 
one consistent with a B~ decay detached from the primary vertex. The invariant mass 
spectrum of D^fi~fi~ candidates is shown in Fig. 10 The mass resolution is 15.5±0.3 
MeV. 
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Figure 10: Invariant mass spectrum for B — > Dffi fi candidates. The line shows the 
fit to the data excluding the B~ mass signal region. 



There are 12 events within the B candidate mass region; it appears that there is 
a dip in the number of events here. An unbinned fit to the data in the sidebands gives 
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an estimate of 22 events. The fluctuation at the B~ mass, therefore, is about two stan- 
dard deviations. Peaking background contributions at the level of current sensitivity are 
negligible (~3xlCT 4 ); thus only combinatorial background is considered. 

After selecting the events in the B~ signal region, we plot the D£fi~ invariant mass 



distribution, which is shown in Fig. 11 A background estimate is made using the sideband 
data in B~ candidate mass (see Fig. 10), by fitting to a 4th order polynomial. The 



background estimated from the sidebands is also shown in the figure. The normalization 
is absolute and in agreement with the data. The data in the signal region is consistent with 
the background estimate. The systematic error due to the fitting procedure is estimated 
using the difference between this fit and the one obtained using a 6th order polynomial. 
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Figure 11: Invariant mass spectrum of -D+/i~ from B~ — > D+{j,~/j,~ events in the signal 
region with the background estimate superimposed (solid curve). There are two combi- 
nations per event. 

The overall efficiencies for different values of are shown in Fig. [8j As done previ- 
ously, during the scan over the accessible Majorana neutrino mass region we use a ±3a^ 
mass window around a given Majorana mass. The resolution is plotted in Fig. [7] as a 
function of Mat. Systematic uncertainties are listed in Table [2] 

Again we provide upper limits as a function of the Majorana neutrino mass, shown in 
Fig. |9^b), only taking into account combinatorial background in this the peaking 

background is absent. For neutrino masses below 5 GeV, the limits on the mass dependent 
branching fractions are mostly < 6 x 10~ 7 . We also determine an upper limit on the total 
branching fraction. Since the background estimate of 22 events exceeds the observed level 
of 12 events we use the CL S method for calculating the upper limit [20]. Assuming a 
phase space decay of the B~ we find 

B[B- D+yrpr) < 5.8 x 1(T 7 at 95% CL. 
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7 Analysis of B 



o„+..-..- 



D U 7T 



A prolific source of neutrinos is semileptonic B~ decay. Majorana neutrinos could be 
produced via semileptonic decays as shown in Fig. |2](c) . Here the mass range probed is 
smaller than in the case of tt + yT /i~ due to the presence of the D° meson in the final state. 
The sensitivity of the search in this channel is also limited by the need to reconstruct the 
D° — > K~n + decay. We do not explicitly veto D* + — > n + D° decays as this would introduce 
an additional systematic uncertainty. The invariant mass distribution of D°7r + /i~ fi~ is 



shown in Fig. 12. The mass resolution is 14.4±0.2 MeV. 
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Figure 12: Invariant mass distribution of D°7i + fi fi 
the data in the sidebands of the B~ signal region. 



The solid line shows a linear fit to 



Peaking backgrounds are essentially absent; the largest source is B — > D°tt it 
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which contributes only 0.13 events in the signal region. The combinatorial background, 
determined by a linear fit to the sidebands of the B~ signal region, predicts 35.9 events, 
while the number observed is 33. 

The 7r + /i~ invariant mass for events within two standard deviations of the B~ mass 
is shown in Fig. [13} The background shape is estimated by a 5th order polynomial fit to 
the sideband data (see Fig. 12) and also shown on the figure. The systematic error on 



this background is estimated using a 7th order polynomial fit. 

The ix + mass resolution is shown in Fig. [7j The dependent efficiencies are 
shown in Fig. [8] They vary from 0.2% to 0.1% over most of the mass range. Systematic 
errors are listed in Table [2] The largest sources of error are the trigger, and the 
dependent efficiencies. 

The upper limits for B{B~ — > D°tt + ^t) as a function of the ir + fi~ mass are shown 
in Fig. [91(c) . For Majorana neutrino masses < 3.0 GeV, the upper limits are less than 
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Figure 13: Invariant mass distribution of 7r + /i~ for B~ — > D°^T ^T-k + in the signal region 
and with estimated background distribution superimposed. There are two combinations 
per event. 



1.6x10 6 at 95% CL. The limit on the branching fraction assuming a phase space decay 



is 



B(B~ D n + /j-fi-) < 1.5 x 1(T 6 at 95% CL. 



8 Conclusions 

A search has been performed for Majorana neutrinos in the B~ decay channels, 
£>(*)+ 

ji jj, , 7t + n jj, , Dfjj, ft , and D°n + fi ji that has only yielded upper limits. The 
ji jj, channels may proceed via virtual Majorana neutrino exchange and thus are 
sensitive to all Majorana neutrino masses. They also could occur via the same annihila- 
tion process as the other modes, though this would be Cabibbo suppressed. The other 
channels provide limits for neutrino masses between 260 and 5000 MeV. The bounds are 
summarized in Table HI These limits are the most restrictive to date. 

Our search has thus far ignored the possibility of a finite neutrino lifetime. Figure [T4| 
shows the relative detection efficiency as a function of Majorana neutrino lifetime, for (a) 
B~ — > [T for a mass of 3 GeV, (b) B~ — > Dffi~fi~ for a mass of 3 GeV, and (c) 
B~ — > D°7r + fi~ fi~ for a mass of 2 GeV. All sensitivity is lost for lifetimes longer than 
10~ 10 s to 10 -11 s, depending on the decay mode. Note that for the ^ final states 
the detection efficiency is independent of the neutrino lifetime, since the neutrino acts a 
virtual particle. 

Our upper limits in the tt + final state can be used to establish neutrino mass 
dependent upper limits on the coupling |V^ 4 | of a heavy Majorana neutrino to a muon and 
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Table 4: Summary of upper limits on branching fractions. Both the limits on the overall 
branching fraction assuming a phase space decay, and the range of limits on the branching 
fraction as a function of Majorana neutrino mass {Mm) are given. All limits are at 95% CL. 

Mode B upper limit Approx. limits as function of 
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Figure 14: Relative efficiencies as a function of Majorana neutrino lifetime for (a) B~ — > 
7r + /i~/i~ for a mass of 3 GeV, (b) B~ — > ' \T for a mass of 3 GeV, and (c) B~ — > 
D°7r + /i~ fi~ for a Majorana neutrino mass of 2 GeV. Where the error bars are not visible, 
they are smaller than radii of the points. 



a virtual W. The matrix element has been calculated in Ref. 0. The results are shown 
in Fig. 15 as a function of M N . A model dependent calculation of B{B~ — > D°rr + /i~ /i~) 
can also be used to extract |V^ 4 | [10J, but the mode is more sensitive. For 

the piT channels upper limits cannot be extracted until there is a theoretical 

calculation of the hadronic form-factor similar to those available for neutrinoless double 
(3 decay. 
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